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K^ ■ Abstract 

^v^ I We use a computer algebra system to compute, in an efficient way, 

[~^ ■ optimal control variational symmetries up to a gauge term. The symme- 

("^ ' tries are then used to obtain families of Noether's first integrals, possibly 

^T , in the presence of nonconservative external forces. As an application, we 

^^ ■ obtain eight independent first integrals for the sub-Riemannian nilpotent 

problem (2, 3, 5, 8). 
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^ ■ 1 Introduction 



The concept of variational symmetry entered into optimal control in the seven- 
ties of the twentieth century [3j. Variational symmetries, which keep an optimal 
control problem invariant, are described mathematically in terms of a group of 
parameter transformations: two transformations performed one after another 
may be replaced by one transformation of the same family; there exists an iden- 
tity transformation; to each transformation there exists an inverse one. Vari- 
ational symmetries are very useful in optimal control, but unfortunately their 
study is not easy, requiring lengthy and cumbersome calculations J15j . 
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Recently there has been an interest in the apphcation of Computer Algebra 
Systems to the study of control systems, and collections of symbolical tools are 
being developed to help on the analysis and solution of complex problems. The 
first computer algebra package for computing the variational symmetries in the 
calculus of variations was given in Q ; then extended to the more general setting 
of optimal control ■ 

In this work we provide a new Maple package for the automatic computation 
of variational symmetries and respective Noether's first integrals in the calculus 
of variations and optimal control. The present package generalize the previous 
results in j2| by introducing two new possibilities: (i) invariance symmetries up 
to a gauge term jjl]; (ii) presence of nonconservative external forces |H|. More- 
over, the efficiency in computing the variational symmetries is largely improved 
when we compare the running times with the ones in [Sj. With the improve- 
ments in the efficiency of the package, we are now able, for the first time in the 
literature, to obtain eight independent first integrals for the nilpotent problem 
(2, 3, 5, 8) of sub-Riemannian geometry. 

2 Nonconservative forces 

Without loss of generality, we consider the optimal control problem in Lagrange 
form: to minimize an integral functional 

/[x(.),u(.)]- / L(i,x(t),u(t))dt (1) 

J a 

subject to a control system described by a system of ordinary differential equa- 
tions of the form 

i(t) = (p(i,x(t),u(i)), (2) 

together with appropriate boundary conditions, not relevant for the present 
study (the results of the paper are valid for arbitrary boundary conditions). The 
Lagrangian L : M x M" x M" ^ M and the velocity vector y. : M x R" x M™ ^ M" 
are assumed to be continuously differentiable functions with respect to all their 
arguments. The controls u: [a, 5] ^ fi C R™ are piecewise continuous functions 
taking values on an open set 17; the state variables x: [a, h] —> R" continuously 
differentiable functions. 

The resolution of optimal control problems usually goes by identifying the 
Pontryagin extremals ^U] . In presence of nonconservative external forces F : R x 
R" X R"* ^ R" the Pontryagin Maximum Principle (PMP) takes the following 
form |S]. 

Theorem 1 (PMP under a nonconservative force F). // (x(), u()) is a 
solution of the optimal control problem (Q)-(0) under presence of a nonconser- 
vative force F (f,x, u), then there exists a non-vanishing pair ('lpo,^p{■)), where 
tpo < is a constant and ip{') o, n-vectorial piecewise C^ -smooth function with 
domain [a,b], in such a way the quadruple (x(-), u(-), i/iq, ■0(')) satisfy the fol- 
lowing conditions almost everywhere in [a, 6] .' 



(i) the nonconservative Hamiltonian system 

fx(t)T = ||(i,x(i),u(i),i^o,V'W), 

\m^ = -^{t, x(t), u{t), V.0, t/'W) + F(i, x(i), u(i))^ 

(ii) the maximality condition 



(3) 



H{t, x(i), u(i), 1^0, ^W) = max i/(i, x(i), v, 1^0, V'W) ; (4) 

where the Hamiltonian H is defined by 

H{t, X, u, iZ-o, ■0) = V'o^(i, X, u) + ■j/)'^ • <p(t, X, u) . (5) 

Remark 2. T/ie right-hand side of the equations of the nonconservative Hamil- 
tonian system ^ represent a row-vector. First equation in ^ is nothing more 
than the control system 0); the second equation is known as the nonconservative 
adjoint system. 

Definition 3. A quadruple {x{-),u{-),tpQ,xp{-)) satisfying Theorem^is said to 
be a nonconservative extremal. A nonconservative extremal is said to be normal 
when ipQ 7^ 0, abnormal when ■(/^o = 0. 

Remark 4. Since we are assuming 51 to be an open set, the maximality condi- 
tion (0J implies the stationary condition 

f)J-f 

— (i,x(t),u(t),Vo,'0(i)) = O, te[a,b]. (6) 

3 Invariance up to a gauge term 

Let h" : [a, 5] X R" X M™ X R X M" ^ M X M" X R™ x M" be a one-parameter 
group of C^ transformations of the fornr 

h''(i,x, u, Vo,'0) = 

(ftj (i, X, u, Vo, lA), hj(i, X, u, Vo, ■0), hu(t, X, u, -0o, -0), h^(t, x, u, ^0, ■0)) • (7) 

Without loss of generality, we assume that the identity transformation of the 
group O is obtained when the parameter s is zero: 

h°{t, X, u, tpo, tp) = t, hl{t, X, u, ipo,tp) = X, 
h°(i,x,u,'(/'o,i/') =u, h^(t,x,u,'(/'o,i/') =■0- 

Associated with a one-parameter group of transformations (|7|l , we introduce its 
infinitesimal generators: 



d 

OS 

d 

U(t,x,u,'(/'o,'</') = -Q-h.1 



d 

, X(i,x,u,V'o,'»/') = -^K 

d 

, *(i,x,u,Vo,'</') = T;-h^ 



(8) 



Definition 5 (Invariance up to a gauge term). An optimal control problem 
Qj-0j is said to be invariant under a one-parameter group of transformations 
Q) up to a gauge term g''{t,yi,u,iPo,ip) e Ci([a, 6],M",M™,R,M";M), if for all 
s sufficiently small and for any subinterval [a, f3] C [a, b] one has 

' H{t^^^{t^),u^t^),^o,r{n) -'^'in^ ~^'in) ^t^ 



H{t, x(i), u(i), 7/^0, t/'W) - rp{tf ■ -x(i) 



+ ^g^(i, x(t), u{t), V'o, tPit))] dt , (9) 

where a'* = /i((a,x(a), u(a),-!/'o, '0(q^)); Z?'* = hf{P,x{P),u{(3),ijJo,tp{P)), and 

When we write Q in terms of the generators ^ , one gets a necessary and 
sufficient condition of invariance - cf. O E] . 

Theorem 6 (Necessary and sufficient condition of invariance). An opti- 
mal control problem is invariant under © up to G{t, x, u, i/iq: V') — x"^'*!^: ^i u, ipo,ip)\ _„ 
or, equivalently, ^ is a symmetry of the problem, up to G, if, and only if, 

ot ox ou dip dt dt dt 

with H the Hamiltonian (|5|). 

Remark 7. The function G(t, x, u, -00, ■0) = g^f;'*(i,x, u,'0Oi''/')| _o *'^ ^^'^'^ 
known in the literature as a gauge term. 

Proof. Transforming the integral on the left-hand side of 1^1 to the interval 
[a,/3], and having in mind that (jHl is satisfied for all subintervals [a,/3] C [a,b], 
the invariance condition can be written in the following equivalent form: 

( dh° (t.x,u,Vio,V)) \ TisU I l\ 

Hih^it,x,u,i.,,iP))-hUt,x,u,^o,^Pf.^^^^^^ 

= H{t, X, u, 00, V) - ■0'^ • ^x + -^g'it, X, u, 00, '0) • 
Differentiating both sides of the equation with respect to s, 



d_ 
ds 



^ dh° (t,x,u,V)o,i/)) 

if (h^(i, X, u, 00, tA)) - h;(i, X, u, 00, V-)^ • ^xg*;!!^^ 

^ dt 



dht{t,x,u,tpo,ip) 
"" dt 



ds ['^9'{t,^,u,fPo,i') ] , 



we obtain the equality 

T dhJ/dA d d/if (dH(\i'')dhl dH{h^) dh^ 



i/(h") - h: 



V 



dhf/dtj dt ds \ dhl ds ah J ds 

dHih") dK dHih") dh'^ dh^^ dhj/di 



9h^ ds ah^ ds ds dhl/dt 



i' I d^ d^dh£ ^f dt ds 

dt dt dt / / 



Finally, choosing s = 0, we express the condition in terms of the infinitesimal 
generators (|HJ| and the fimction ^^(ijX, u, ■00,1/') = ■^g''{t,x,u,1|JQ,^p)\ _ : 



D 



dt V dt ax 



9i/ ^, OH ^ 


*'^x 


ou dip 





,T / dX .dT\\ dG 



4 Nonconservative Noether's theorem 

Emmy Noether was the first who established the relation between the existence 
of invariance transformations of the problems and the existence of conservation 
laws - first integrals of the Euler-Lagrange or Hamiltonian equations [S] . A gen- 
eralization of the classical result of E. Noether for the nonconservative calculus 
of variations was recently given by Fu and Chen 6 ; then extended to the more 
general setting of optimal control by Frederico and Torres [5] . 

Using (PJ, together with the stationary condition ©, one can deduce that 
along the nonconservative Pontryagin extremals (Definition O, the total deriva- 
tive of the Hamiltonian with respect to the independent variable t is equal to 
its partial derivative plus the scalar product of the velocity vector with the 
resultant nonconservative forces F [S]: 

-H{t, x(t), u(t), ^oMt)) = -QiH{t, x{t),u{t),yjo,tp{t))+i{tf-F{t, x(t), u{t)) . 

(11) 
Using this fact, the nonconservative optimal control version of E. Noether's 
theorem is easily obtained from the necessary and sufficient invariance condition 
ifTUIl . restricting attention to the quadruples (x(-), u(-), -00, ■»/'(•)) that satisfy 
the nonconservative Hamiltonian system and the maximality condition Q: 
along the extremals, equalities ((SJ, ®, and Ullfl permit to simplify (|fO(l to the 



form 



dt J \ ^ J ^ dt dt dt 

dH dT -T T dX dG , T tn 

^ -—T + H—- -tp -X-ip^ ■- —- (±^T- X^) • F 



^-^(i/T-^'^-X-G- /(i'^T-X^) -Fdi) =0. 



dt dt dt dt 

d 

dt 

This means that HT - ip'^ ■ X - G - j{^^ T - X^) ■ Fdt is a first integral 
whenever the optimal control problem under consideration admits a symmetry 
^ up to the gauge term G: 

Theorem 8 (Nonconservative Optimal Control version of Noether's 
Principle). // the infinitesimal generators ^ constitute a symmetry of the 
optimal control problem OP-0' under presence of nonconservative forces with 
resultant vector F (t, x, u), then 

±{tfTit, x(t), u(t), Vo, t/^W) -X(t, x(t), u(t), ^o,rPit)f) -Fit, x(t), u(t))dt 

+ iPitf ■ X(t, x(t), u(t), iJo,tp{t)) + G{t, x(t), u(t), 7/;o, -0(t)) 
-H{t, x(t), u(t), 1^0, ^{t)) T{t, x(t), u(t), 7/^0, '0(t)) = const 

(12) 

is a conservation law, i.e., condition (|12|l holds for all t in [a,b] and for every 
nonconservative extremal (x(-), u(-), -00, '!/'(■)) of the problem. 

5 Computation of symmetries up to a gauge term 

The main problem in obtaining Noether's conservation laws (in applying The- 
orem |H} resides in the determination of the symmetries and respective gauge 
terms. If n effective first integrals exist ^21: then the optimal control prob- 
lem is integrable, and classical results allow the integration of the equations of 
motion. 

Here we propose an algorithm for determining the infinitesimal generators |(SJ) 
and the gauge terms G which define a variational symmetry. Let us assume, for 
the moment, that the optimal controls are G^ functions (in iJ7|we will drop this 
restrictive assumption, just by assuming that T, X, and G do not depend on the 
control variables). The key point to compute symmetries consists in generalizing 
the method used in |S1 §3] to the nonconservative and gauge-invariant cases. The 
idea is simple: when we substitute the Hamiltonian H and its partial derivatives 
in the invariance identity H1(J|I . then the condition becomes a polynomial in x, 
li and ip, and one can equal the coefficients of the polynomial to zero. Thus, 
given an optimal control problem Q-ldl, defined by a Lagrangian L and a 
velocity vector ip, we determine the infinitesimal generators T, X, U and ^ and 



the gauge term G, which define a symmetry for the problem, by the following 
method: (i) we define the respective Hamiltonian 10); (ii) we substitute H and 
its partial derivatives into l(Tn|l : (iii) expanding the total derivatives 



(13) 



we write equation (|10|) as a polynomial 

A{t, X, u, Vo, V')+5(^, X, u, i/jQ, 4>)-x+C{t, X, u, '00, ■0)-u+L'(i, x, u, V'o, ■0)-»/' = 

(14) 
in the 2n + m derivatives x, u and ip: 



dT 
dt 


dT dT 
dt ^ dx 


. dT 
du 


. dT ■ 

u+ —- -tp, 
dtp 


dX 
dt 


dX dX 
dt ' dx 


. dX 
9u 


. dX ■ 

■u+ —- -tp 
dtp 


dG 
dt 


dG dG 
dt '^ dx 


. dG 
du 


. dG ■ 

U+ —J -tp, 
dtp 



fdH dH 
\dt dx 

dT 



X 



dH 



U 



dH 



* 



H- 



dy 



tP 



9u 9'0 

T ax 5G' 



* + iJ 



dT 
'dt 



tp^ 



9> 



d 



X 



X 



H 



dT 



H—-tP 
dtp 



du 

T dX 



tP 



dX dG 

~dt~~dt 

T ax dG 



du 
dG_ 

9'0 dtp 



du 

-0 = 



(15) 



The terms in (|15|l . which involve derivatives with respect to vectors, are ex- 
panded in row-vectors or in matrices, depending, respectively, if the function is 
a scalar or a vectorial one. For example. 



dT 


' dT dT 


dT' 










dx 


dxi 9x2 


dXn_ 


; 












r dx^ 


dXi 


dXi 








diPi 


dip2 


dip„ 


dx 


' dX dX 


dx 




dX2 


dX2 

dTp2 


dX2 


dtp 


.5-01 9-02 


dlpn 
















dX„ 


dX„ . 


. ax„ 










L dipi 


dip2 


dip„ 



Equation H15|) is a differential equation in the 2n + m + 2 unknown functions T, 
Xi, . . . , Xn, Ui, . . . , Um, 5^1, ■ ■ • , 'Pn and G. This equation must hold for all 
ii, . . . , in, ui, • ■ ■ , Un, 0'i, ■ • ■ , 0'n, and therefore the coefficients A, B, G and 



D of polynomial p4|) must vanish, that is, 

dt dx du dxp dt dt dt 

ox ox ox 

dT ,T 9X dG 
Ou ou Ou 

ar T ax dG 

dxj} dij} dip 

(16) 
System of equations (|16|l . obtained from p5|) . is a system of 2n + to + 1 partial 
differential equations with 2n + to + 2 unknown functions so, in general, there 
exists not a unique symmetry but a family of such symmetries. The system (|16|l 
becomes even more under-determined when one assumes, as in [jT] that T, X, 
and G do not depend on the control variables u. Although a system of partial 
differential equations, solving H16|) is possible because the system is of the first 
order and linear with respect to the unknown functions and their derivatives. We 
solve the system of PDEs by the method of (additive) separation of variables, 
as explained in j2]. Following '5', the generators are replaced by the sum of 
unknown functions, one for each variable. For example, T(t,Xi,X2,'>pi,'>p2) = 
Tiit) + T2ixi) + T3(x2) + T^iipi) + Tz{^2)- When dealing with optimal control 
problems with several state and control variables, the number of calculations is 
big enough, and the help of the computer is more than welcome. We define a 
Maple procedure Symmetry that does all the cumbersome calculations for us. 
The procedure receives, as input, the Lagrangian and the velocity vector; and 
returns, as output, a family of symmetries (T, X, U, >&) and, if necessary, the 
respective gauge term G. We remark that since system (|16|l is homogeneous, 
we always have, as trivial solution, (T, X, U, >&) = 0. 

6 The computer algebra package 

We obtain Noether conservation laws, in an automatic way, through two steps: 
(i) with our procedure Symmetry we obtain the variational symmetries and re- 
spective gauge terms; (ii) using the obtained symmetries, gauge terms, and non- 
conservative forces as input to procedure Noether^ we obtain the correspondent 
conservation laws. In ^we give several examples, not covered by the previous 
results in UllS], illustrating the whole process. Given the limit on the maxi- 
mum number of pages of the paper, we do not provide the Maple definitions for 
the procedures Symmetry and Noether here. The complete Maple package can 
be freely obtained from http://www.mat.ua.pt/delfim/maple.htm together 
with an online help database for the Maple system. 

Novelties of the procedures Symmetry and Noether with respect to the pre- 
vious versions in |7||H| are: (i) possibility of procedure Symmetry to cover invari- 
ance symmetries up to a gauge term, according with iJ3£^iid i^l (ii) improvements 



of efficiency - see [El (iii) possibility of procedure Noether to consider problems 
of the calculus of variations and optimal control under nonconservative external 
forces, according with 21 (iv) improvement of the usage of the procedures by 
introduction of several optional parameters, as illustrated in ^ Moreover, a 
new Maple procedure called PMP was added which implements Theorem ^ ac- 
cording with 32^ The procedure PMP is very useful in practice, when dealing 
with concrete problems of the calculus of variations and optimal control - cf. 
^ The input to the procedure is: the Lagrangian L and the velocity vector cp, 
that define the optimal control problem iQJ-© and the respective Hamiltonian 
H; the nonconservative external forces (if present); and several useful optional 
arguments which define the output. The output of PMP is either (depending 
on the optional parameters): the (nonconservative) extremals; the equations of 
the (nonconservative) Hamiltonian system and stationary condition; or, alter- 
natively, the Hamiltonian. We refer the reader to the Examples on iJHl for a 
general overview on the usage of the developed Maple procedures; to the an- 
notated Maple worksheet available at http://www.mat.ua.pt/delfim/maple.htm, 
with all the definitions of the package, detailed documentation, and many other 
examples not given here, for more details. The reader is free to experiment 
the Maple package in order to determine variational symmetries and Noether 
conservation laws on his/her own problems. 

7 Efficiency, comparison with previous results 

The high number of dependences that the infinitesimal generators may present, 
affect, excessively, the efficiency of the method described in ^ namely for prob- 
lems with a large number of state and control variables. In order to quantify 
this effect, we measured the computing running times of our procedure Symme- 
try for different dependences of the infinitesimal generators lO, with a large set 
of optimal control problems: the ten problems considered in [HI §4, §5] (exam- 
ples 4.1-4.6 and 5.1-5.4), together with twelve new problems. Three of these 
new problems are given in ^Hl the complete set of problems being available as 
a Maple worksheet, as mentioned in ^ All the computational processing was 
carried out with the Maple 10 Computer Algebra System on a 1.4GHz Pentium 
Centrino with 512MB of RAM. In the previous work |21, the maximum number 
of dependences for each generator, as indicated in ©, is always considered. We 
denote here such situation by Dl. In the Dl case, and as noticed in |S], the 
involved computational effort is sometimes very high: the computing times in- 
crease exponentially with the dimension of the problem. This is particularly well 
illustrated with the following problems of sub-Riemannian geometry: nilpotent 
problem (2, 3), with three state variables, requires a total computing time of one 
minute ([HI Example 4.5]); problem (2,3,5), with five state variables, requires 
thirty minutes (jSl Example 4.6]); the problem (2,3,5,8), with eight state vari- 

^In the software Cotcot, available from http://www.n7.fr/apo/cotcot/, the tool Adifor for 
automatic differentiation of Fortran is also used to generate, in the conservative case, the 
equations of the Pontryagin maximum principle 0. 



ables, was not studied in [5], and thought to be out of its capacities. We compute 
here its symmetries in Example II II with the present Maple package, with forty 
one minutes of computing time; while the method in [5] requires, approximately, 
thirty times this value: twenty hours of computing time are needed.^ 

The computing running times largely depend on the numbers n and to, re- 
spectively the number of state and control variables: besides directly influencing 
the number of dependences of the unknown functions (infinitesimal generators), 
they determine the amount of those functions and the number of partial differ- 
ential equations that must be solved in order to find the variational symmetries. 
Without considering the gauge term, we come across a system oi m + 2n + 1 
partial differential equations and m + 2n+l unknown functions, each one of the 
unknown functions being dependent of m -|- 2n + 1 variables. We address here 
the following question: is there some way to simplify the process of obtaining 
the variational symmetries? 

Although knowing that the complexity of the method is intimately related 
with the values n and to, that are fixed with a given optimal control problem, 
we get, even so, a quite satisfactory answer to the question. Analyzing the 
results from the test set of problems, we verify that, in spite of considering the 
maximum number of dependences {Dl), the infinitesimal generators obtained 
through the procedure Symmetry are, nevertheless, almost always, dependent 
functions of a quite reduced number of variables. When we restrict ourselves 
to the dependences T{t), X(i,x), \J{u,^p), '^{^p) - that we identify as D2 - 
we are able to cover the totality of the twenty two considered problems in our 
study. If in the formulation of the system of PDEs H16(l we only enter with these 
dependences, besides the obvious reduction of the number of dependences of the 
unknown functions, we reduce the number of equations for less of half: from 
771 + 2n-|- 1 to n+ 1. In agreement with the simulations done, the efficiency of the 
procedure Sym,m,etry increases significantly with this new group of dependences 
{D2). For instance, for the problem (2,3,5) of sub-Riemannian geometry |S1 
Example 4.6], a problem with two controls and five state variables, the running 
time passed from half an hour to less than one and a half minute. We have also 
considered another more simplified set of dependences, denoted by D3: T{t), 
X(i,x), U(t, u), '4'(t, -0). With it, is now possible to obtain the symmetries of 
the sub-Riemannian nilpotent problem (2, 3, 5, 8) f Example 1 11 |l . in less than 45 
minutes; and it is still possible to obtain the same conservation laws for all the 
twenty two studied problems (in three of the problems |H1 Examples 4.4, 5.2 and 
5.3] the generators were different, since the more general generators U depend 
on the variables ip, but the correspondent Noether conservation laws H12|) are 
exactly the same since they only depend on the generators T and X). Finally, 
we repeated the study for a more restricted group of dependences (-D^): T{t), 
X(x), U(u), ^t'(■0). As expected, the time of processing suffered an additional 
reduction (for the (2, 3, 5, 8) problem the running time passed from 44'16" to 
28'21"), but in this case not all the family of conservation laws for the problems 

^We believe that the forty minutes of computing time can still be diminished by us- 
ing a programming language closer to machine, for instance using Adifor: http://www- 
unix.mcs.anl.gov/autodiff/ADIFOR. 
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are obtained. For four of the problems - [S] Example 4.3], Examples 2 and 3 
in the Maple worksheet, and Example]^- only particular cases of the complete 
family of conservation laws are obtained. 

To summarize the influence that the different dependences of the generators 
have in the efficiency of the procedure Symmetry, we give in Table ^ the run- 
ning times for computing the variational symmetries of the three problems of 
sub-Riemannian geometry already mentioned: [3 Examples 4.5 and 4.6] and 
Example E] All the three problems have two control variables and the same 
Lagrangian, but a different number of state variables, respectively, 3,5, and 8. 



Dependences 


N° PDEs* 


prob. (2, 3) 


prob. (2, 3, 5) 


prob. (2, 3, 5, 8) 


Dl 


m+2n + l 


1'04" 


30' 34" 


20/i07'12" 


D2 


n + 1 


5" 


1'26" 


51'28" 


D3 


n + 1 


4" 


1'09" 


44' 16" 


D4 


n + 1 


2" 


38" 


28'21" 



* n = n° of state variables; m = n° of control variables. 

Table 1: Running times of procedure Sym^metry for three problems of sub- 
Riemannian geometry (|Si Examples 4.5, 4.6] and Example 111(1 . with differ- 
ent dependences of the infinitesimal generators: Dl - [T(t,x,u, ■!/>), X(i,x,u, ■;/)), 
U(f,x,u,^), *(f,x,u,V)]; D2 - [T(i), X(t,x), U(u,^), *(V')]; D3 - [T{t), X(i,x), 
U(t,u), *(t,V)]; D4 - [T{t), X(x), U(u), *(V')]. 



We verify that of the four sets of studied generators, just with D4 it was not 
possible to obtain, with full generality, the totality of Noethcr's conservation 
laws for the twenty two considered problems. The set of generators D3 {T{t), 
X(i,x), U(i, u), \I/(t, ■0)) gives the best compromise: it presents the best run- 
ning times, between the generators that give the complete family of variational 
symmetries and Noether conservation laws for the problems we have studied; 
running times arc much better than the ones obtained with the generators Dl. 
We recommend the user to try configuration D3 first on his/her own optimal 
control problems. Considering t and x for the dependences of the gauge term 
- G{t, x) - the system of PDEs that we have to solve, in order to find the 
variational symmetries, takes form (cf. H16|l ) 



'di 



^' 



■T 



dH 



X 



dH 

du 



U- 



dH 



^ + H 



'dt 



^P' 



'dt 



dG 

'at 



= 0, 



ax dG 



ax ox 



(17) 

Our present procedure Symmetry computes, by default, the variational symme- 
tries as defined by D3, and with a gauge term G{t,x): by default Symmetry 
solves system (|17|l . Through optional parameters, it is possible to find the 
variational symmetries for other generators and gauge terms: in order to use 
all the dependences {Dl ) one must use option alldep; to use a minimum of 
dependences {D4-) one uses option mindep. We remark that with the class of 
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generators D3, T and X are not functions of u, and there is no need to assume 
the control variables u to be smooth functions (cf. (|13() ). 

Table 12 shows the computing running times needed to obtain all the varia- 
tional symmetries of the problems in O §4, §5], by using the default version of 
procedure Symmetry we give here (generators D3)] and by using the version in 
|S], which is a particular case of our present procedure - see 10 for examples not 
covered by the previous methods in jHj - obtained using option alldep, that is, 
generators Dl. The time needed to compute the variational symmetries for the 
(2, 3, 5) problem (Example 4.6 in [5]) decreased from thirty minutes to one. 





4.1 


4.2 


4.3 


4.4 


4.5 


4.6 


5.1 


5.2 


5.3 


5.4 


Dl 


2" 


1'13" 


2'44" 


6'41" 


1'04" 


30' 34" 


8" 


17" 


6'42" 


1" 


D3 


0" 


5" 


11" 


18" 


4" 


1'09" 


0" 


3" 


16" 


0" 



Table 2: Running times of procedure Symmetry for all the problems of previous 
work |H], with the generator sets Dl (the only possibility in 0) - [T(t,x,u, ■(/)), 
X(i,x,u,V), U(i,x,u,V), *(t,x,u,i/')], and D3 - [T{t), X(i,x), U(f,u), *(t,V)]. 



The use of generators with a smaller number of dependences leads to a 
drastic reduction of the computing running times. For the studied problems, 
the use of generators D3 permit to obtain the same results while decrease the 
total processing times for about 4% of the ones verified in [5] (generators Dl ) . 

8 Examples of the new possibilities 

In order to show the functionality and the use of the new procedures, we apply 
our Maple package to three concrete optimal control problems which are not 
covered by the previous results in |Z| |Hj. All the examples were carried out 
with Maple version 10 on a 1.4GHz 512MB RAM Pentium Centrino. The run- 
ning time of procedure Symmetry is indicated, for each example, in the format 
min'sec". All the other Maple commands run instantaneously. 

8.1 Variational symmetries up to a gauge term 

We begin with a very simple example of the classical calculus of variations. We 
recall that for the fundamental problem of the calculus of variations there are 
no abnormal extremals, so one can choose ^q = —\ (we use option noabn of our 
Maple package). 

Example 9. (O'OO") Let us consider the following scalar problem of the calculus 
of variations (n = m = \): 



{u{t)) dt — > min , 
x{t) ^u{t). 
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In this case L — v? and tp = u. First we obtain the variational symmetries of the 
problem (Maple procedure Symmetry^ up to a gauge term (parameter gauge). 

> S := Symmetry (u~ 2, u, t ,x,u, showt, gauge ) ; 

r = 2 Cat + Ce, X = i ^ + C2x{t) + C4, t/ = i ^ - u{t)C2, 

* = -V'(t)C2 - C3, GAUGE = C3x{t) + Cs 

Noether conservation laws are obtained through Theorem (Maple procedure 
Noether j with the generators and the gauge term just obtained. 

> CL := Noether (u"2,u,t ,x,u,S,showt ,noabn,H) ; 

CL := (-^ Cst + C2x{t) + CA iP{t) - H (2 Cat + Ce) + C3x(t) + C5 = const 

The Hamiltonian H , which appears in the above family of conservation laws, is 
given by ^: 

> H := PMP(u~2,u,t ,x,u, evalH.showt ,noabn) ; 

H ~ ~u{tf + u{t)iP{t) 

This is a very simple problem, just used to illustrate, in the simplest possible 
way, our Maple procedures. In this case it is an easy exercise to obtain the 
extremals by direct application of the Pontryagin Maximum Principle or the 
Euler- Lagrange equations, 

> extremals := PMP(u~2,u,t ,x,u,showt ,noabn) ; 

r 1 1 

extremals := < tp{t) = K2, x{t) = - K2 t + A'l, u{t) = - K2 

and one can validate the obtained conservation laws by applying the definition of 
conservation law: by definition, the obtained family of conservation laws must 
hold along all the extremals of the problem. 

> subs (extremals, CL) ; 

K2 CiKr + A-a C4 - i KICq + C3K1 + C5 = const 

8.2 Presence of nonconservative forces 

We consider now a problem of the calculus of variations under the action of a 
nonconservative force. The problem is borrowed from I", §4]. 

Example 10. (ji = 1, in = 2, O'Ol") The problem is defined by the Lagrangian 

L{q,q,q) = h^i^) + h'^'iit) + h^git) , and presence of the nonconservative 

2 
force f{t) — fi q(t) + ^q{t) — 2 ^ Q (t) which depends on higher-order derivatives 

(a, b, and 11 are constants). 

> PDEtools [declare] (prime=t) ; 
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derivatives with respect to t of functions of one variable will now be displayed with 

> L := u"2/2+a*v"2/2+b*q"2/2; 

> phi : = [v , u] ; 

> f := mu*v+inu"2/a~2*u-2*mu/a*z(t) ; 

r -'^2,1 2,172 

■" 2 " 2 "" 2 ^ 

ip ■- [v, u] 

fl^U fizjt) 

f ■- ^uH 5- - 2 

a^ a 

> S := SymmetryCL, phi, t, [q,v] , u) ; 

S --[T^ Ci, Xi=0, X2=0, U = 0, *i = 0, *2 = 0] 

> CL := Noether(L, phi, t, [q,v] , u, S, ncf=[f,0], noabn) ; 

CL := - (-^uitf - ^av(tf - i bq{tf + ^i(i)t;(i) + MtMt)) Ci 



+ |c,g'(^M.(i) + ^-2^ 



dt = const 



The multipliers ipiit) and ip2it) are obtained using the adjoint system and the 
stationary condition, as given by Theorem^ 

> sys := PMP(L, phi, t, [q,v] , u, noabn, evalSyst, ncf=[f,0], showt) ; 

{q' = vit),v' = u{t)} , {-^1' = ^f.v{t) - ^^ + 2 ^ - bqit), 
--02' = -av{t) + ^i(t) \,{-u{t) + ip2{t) = 0} 



sys 



> dsolve({sys[2] [2],sys[3] [] },{psi [1] (t) ,psi [2] (t)}) ; 

{M't)^<'t), i^iit) = ~u +av{t)} 
With substitutions 

> subsC/., z(t)=diff (u(t),t), u(t)=diff (v(t),t), v(t)=dif f (q(t) ,t) , 

C[l]=l, CL); 

1 "2 ,1 '2,1, u\'2 I '" , i\ I , I I ( ' , M^e" r,i^g"'\ ,, , 

"2? +2"^ +2^<iV')-[-<l +aq)q+ q [1^1 +—2 2 —— j dt = const 

one obtains the conservation law 01 §^y. We remark that the conclusion is 
nontrivial, and difficult to obtain without Noether's principle. 
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8.3 The sub-Riemannian nilpotent case (2, 3, 5, 8) 

We finish the section by applying our Maple package to one important problem: 
the study of sub-Riemannian geodesies. The reader, interested in the study of 
symmetries of flat distributions of sub-Riemannian geometry, is referred to |13) . 
Here we use a formulation of the nilpotent problem (2, 3, 5, 8) which is obtained 
using the results of [TT] . 

Example 11. (4-4'16") The problem can be defined in the following way: 

xi{t) = Ul{t) , 
i2it) = U2{t) , 
X?,{t) =U2{t)xi{t), 



(uiitf + U2itf^ dt 



iiit) = ^U2{t)xi{t) , 

ibit) = U2it)xi{t)x2{t) , 

X(i{t) = ^U2it)xi{tf , 

Xj{t) = lu2it)xi{t)'^X2{t), 

^Xs{t) = ^U2it)xi{t)x2{tf ■ 

The integrability of the problem is still an open question Jfl'A \iyi/ . but eight in- 
dependent conservation laws can be determined with our present Maple package. 

> L := l/2*(u[l]-2+u[2]-2); 

>phi: = [u[l], u[2], u[2]*x[l], (u[2] /2)*x[l] "2, u[2] *x[l] *x[2] , 

(u[2]/6)*x[l]-3, (u[2]/2)*x[l]-2*x[2], (u[2] /2)*x[l] *x[2] "2] ; 

> XX := [x[i]$i=l. .8] ; 

> UU := [u[l] ,u[2]] ; 

T 1 2,1 2 

L .= -ui +-U, 



V 



I2 I3I2I 2 

Ul, U2, U2Xi,-U2Xl , U2XlX2,-U2Xl , - U2X1 X2,-U2XlX2 
2 2 A 



XX ~ [xi,X2,X3,X4,X5,X6,Xy,Xs] 

UU :— [ui, U2] 
> SymmetryCL, phi, t, XX, UU) ; 



T — Cit + Ct, Xi — - Cixi, X2 = C2 + - C1X2, X3 = C1X3 + Cg, 

3 3 

X4 = - C1X4, + C'e, Xs — C2X3 + - Cixs + C3, Xq — 2 C'lXa + C5, 



Xr = C2X4 + 2 CiX7 + Cg, Xs = C2X5 + 2 Cixs + C4, Ui 



1 n 



U2 



■C1U2, *i 



1 



ClV'l,*2 



■ Cl^2, *3 = -1P3C1 - C2ll)5, 



*4 



■V4C1 -C2-tp7, *5 



■C\^5-C2^8, *6 = -2ClV'6, 



*7 



-2ClV'7, *8 = -2ClV8 
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> CL := NoetherCL, phi, t, XX, UU, "/., H) ; 

CL ■- - CixiTpx + f C2 + - C1X2 \ip2 + {C1X3 + Cs) V's + ( 2 C!iXA + Ce j i/'4 

+ [ C2X3 + - C1X5 + C3 j Vs + (2 C1X6 + Cs) V'6 + (C2X4 + 2 Cia;7 + C9) i>7 
+ (C2a;5 + 2 Cia;8 + C4) ips - H {Cit + C7) = const 
The Hamiltonian is given by 

> Hamilt := PMP(L, phi, t, XX, UU, noabn, evalH) ; 



Hamilt ■— --Ui^ - - U2^ + IplUi + tp2U2 + i>zU2Xx + - '4'aU2X\^ + iIJ5U2XxX2 

+ - U2Xl^^l)e + - U2Xx^X24>7 + - U2XlX2^i)S 
O 2 2 

and the extremal controls are obtained through the stationary condition. 

> PMP(L,phi,t, XX, UU, noabn, evalSyst) [3] ; 

-U2 + ■02 + '(/'sa;! + - ipAXi^ + ^p5XiX2 + - xi^^pe + - a::i^a;2'!/'7 + - xyX2^'4>s = 0, 

> solve (y.,-[u[l] ,u[2]}-); 

Ui = ^pi, U2 = '4>5XlX2 + ■!/'2 + i>3Xx + - IpAXl^ + - Xl'xpe, + - Xl^ X2i)7 + - X\X2^ IpH 

> H = expaiid(subs(°/o, Hamilt)); 

-W = - ■!/'2a;ia;2^-!/'8 + iphX\X2i>2 + ^5a:i^a;2V'3 + - i>2ipAXi^ + - tpzXi'ipA 

,1,2 2 2,1, 3,,1 2 4,2,1 4 2,2,1, 5 

+ -W5 XX X2 + -W2X\ ¥'6 + -a;i X2 ip» + ^ ^1 ^2 07 + T^r V'4a;i 06 

z Do o iz 

+ - ip3^xi^ + — a;i'^V6^ + o ^2^ + - i/"!^ + - V'4^a;i'' + - i/'5a;i''a;2i/'6 + 7: ■!/'3a;i^a;2^-!/'8 
z /z z z o D z 

+ -TpAXl X2 ■!/'8 + T^4a;i X2'll)7 + ^2^3Xl + -X\ X2 V'tV'S + T^ a^l 1pfiX2'lp7 

+ — x\'^iljQX2^i)a + - ^2a;i^a;2^7 + - ip^xi^ X2ipA + - ■!/'5a;i^a;2^'!/'8 + ^ ■!/'5a;i^a;2^'!/'7 

+ ^5 '4>iXt'X2i>7 + - IpzXl'^XpQ 

Now, the eight conservation laws, we are looking for, are easily obtained: 
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> subs (C [8] = 


1, 


seq(C[i]=0,i=l. .9), 


CD 




> subs (C [6] = 


1, 


seq(C[i]=0,i=l. .9), 


CD 




> subs (C [3] = 


1, 


seq(C[i]=0,i=l. .9), 


CD 




> subs (C [5] = 


1, 


seq(C[i]=0,i=l. .9), 


CD 




> subs (C [9] = 


1, 


seq(C[i]=0,i=l. .9), 


CD 




> subs (C [4] = 


1, 


seq(C[i]=0,i=l. .9), 


CD 




> subs (C [2] = 


1, 


seq(C[i]=0,i=l. .9), 


CD 




> subs (C [7] = 


-1, 


seq(C[i]=0,i=l. .9) , 


CD 








ip3 = const 








tp4, = const 








tp5 — const 








tpe — const 








tp7 — const 








tps — const 








tp2 + Xstps, + 


X4'>p7 + xs'ips = const 






H — const 







Given the results of jll| , one can say that the sub-Riemannian nilpotent Lie 
group of type (2,3,5,8) has seven trivial first integrals: the Hamiltonian H; 
and the multipliers -03, V'4j "051 "06, ^7, '08 ■ Together with the non-trivial first 
integral 02 + 2^305 + x^-ipr + a^508, here first obtained, it is possible to prove that 
the system is integrable. This is nontrivial since Liouville theorem does not 
apply: the set of first integrals is not involutive (for instance, Poisson bracket 
between 03 and 02 + x^ip^ + x^ipj + x^ips is not zero). This question is under 
study and will be addressed in a forthcoming publication. 
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